Effect of the Kondo correlation on thermopower in a Quantum Dot 
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In this paper we study the thermopower of a quantum dot connected to two leads in the presence of 
Kondo correlation by employing a modified second-order perturbation scheme at nonequilibrium. A 
, simple scheme, Ng's ansatz [Phys. Rev. Lett. 76, 487 (1996)], is adopted to calculate nonequilibrium 

■ distribution Green's function and its validity is further checked with regard to the Onsager relation. 

^Sj ' Numerical results demonstrate that the sign of the thermopower can be changed by tuning the 

energy level of the quantum dot, leading to a oscillatory behavior with a suppressed magnitude 
O ' due to the Kondo effect. We also calculate the thermal conductance of the system, and find that 

\ the Wiedemann-Franz law is obeyed at low temperature but violated with increasing temperature, 

^^^^^ . corresponding to emerging and quenching of the Kondo effect. 
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PACS numbers; 72.10.Fk, 72.15.Qm, 72. 20. Pa, 73.23.Hk 



I. INTRODUCTION 



Observation of the Kondo effect in a semiconductor quantum dot (QD)^'^ which provided a testing ground of the 
^ ^ quantum behavior of electron wave functions and many-body effects, has stimulated a great deal of experimental'^'^ 
d , and theoretical^"^ investigation. So far, many interesting features in QD, such as Kondo-assisted enhancement 
^ ■ of conductance, its specific temperature dependence, a peak splitting in a magnetic field, zero-bias maximum of 
I differential conductance in the Kondo regime, and singlet-triplet Kondo effect have been successfully explored in 
' O ■ early theoretical works, which described the electron transport through QD using the well-known impurity Anderson 
' model. 

A convenient way of describing transport through an Anderson impurity under out-of-equilibrium condition is to 
employ the nonequilibrium Green's function (GF) technique of Keldysh. In the framework of nonequilibrium GF, 
one needs a consistent calculations of retarded (advanced) and distribution (lesser) Greens's functions in order to 
study the out-of-equilibrium properties, while only the retarded GF is needed when studying equilibrium properties 
of the system. In literature, the single-impurity Anderson model has been extensively studied for more than thirty 
years and many numerical and analytical methods have been developed to explore its equilibrium properties, such 
^vq ' as equation-of-motion (EOM) combined with the decoupling approximation,^^ the second-order perturbation theory 
. (SOFT) for on-site Coulomb interaction [/,^^ slave-boson non-crossing approximation (NCA).^'^'^'' However, it is 
' still a great challenge to evaluate the distribution GF of the strongly correlated systems under out-of-equilibrium 
04 , condition, because one has to deal with both the on-site Coulomb interaction and the tunnelings between the dot 
and reservoirs simultaneously. To our knowledge, there were three attempts to surmount this technical difficulty. 
The first one has been carried out by Meir,^ who generalized the NCA to the nonequilibrium situation basing on the 
Coleman's slave-boson method in the limit of infinite Coulomb interaction U — > oo.^'^ This scheme transforms the 
strong correlation Hamiltonian into a noninteracting equivalent one by introducing several auxiliary boson operators. 
. Therefore, it evades the dilemma how both the Coulomb interaction U and the mixing between QD and the two 
■ leads can be treated simultaneously in deriving the retarded and lesser (greater) self-energies of QD. It is well known 
O . that NCA provides a good description for the investigation of the excitation spectra of QD, but the necessary Fermi 
O liquid behavior was not reproduced at low energy and low temperature limit. To extend the slave-boson method to 
J> . study Kondo-type transport through double QDs at low temperature, espatially at zero temperature, a slave-boson 
JJ"! ' mean field approach has been presented by replacing the slave-boson operators with their expectation values. Very 
rN ' recently, a finite-fj slave-boson mean field scheme has been developed by us to explore Kondo-type transport through 
> QD^^ and double QDs^^. Nevertheless, these slave-boson mean field approaches are at the same level of approximation 
" " ' in deriving the interacting retarded and lesser self-energies. 

The third attempt has been made by Ng,^ who has developed an approximative scheme, i.e., Ng's ansatz, to 
obtain the lesser GF from EOM. It is impossible to obtain straightforwardly the distribution GF from EOM without 
introducing additional assumptions. Ng assumed that the interacting greater or lesser self-energy of QD is related to 
noninteracting one, and the ratio parameter is determined by the interacting retarded self-energies derived by means 
of EOM. This simple and effective ansatz has later been applied in more complicated structures involved QD, such as 
normal metal-QD-superconductor^*, and superconductor-QD-superconductor^^. Nevertheless, except for these three 
advantages pointed out by Ng in Ref. 9, up to now, no further inspectation about the validity of the ansatz has been 
presented. 
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In this paper wc will provide another validity check of Ng's ansatz by analyzing the particle current and thermal flux 
through QD connecting to two reservoirs on the basis of Anderson single-impurity model by means of noncquilibrium 
GF with the help of Langreth continuation rules and Ng's ansatz. We find that the resulting particle current and 
heat flux driven by bias voltages and temperature gradients between two leads satisfy the Onsager relation in nearly 
equilibrium regime, thus provide a natural check on the validity of derivation of the lesser GF. 

To date, most theoretical calculation and experimental measurement on the thermopower of QD have focused in 
the Coulomb blockade (CB) regime, which have revealed that when the gate voltage defining QD is swept, the 
thermopower oscillates about zero (sawtooth behavior) with a period equal to that of the CB oscillations in electric 
conductance. However, thermopower across QD in the Kondo regime is still much less studied. ^^'^'^ The second 
purpose of this paper is to investigate the Kondo effect on the linear thermopower of QD basing on Ng's ansatz. 

We organize the rest parts of the paper as following. In the second section we derive particle current and thermal 
flux formulas through interacting QD within the framework of noncquilibrium GF. The electric and thermoelectric 
transport coefficients are derived in the presence of both chemical potential and temperature gradients between two 
leads and automatically satisfy Onsager relation in the linear transport regime. We note that within this approx- 
imation, the same current formula as that in Ref. 25 is derived without the presumption of proportional coupling 
ri(w) = XTii{Lo). In section 3, numerical calculation of linear Kondo-type thermopower S in QD are reported as a 
function of gate voltage and temperature, which shows a gate- voltage-controlled change of sign and largely suppressed 
oscillatory magnitude due to the Kondo effect. We also discuss the thermal conduction coefficient k and its violations 
of the well-known Wiedemann-Pranz law. Finally, a conclusion is given in Section 4. 



II. THERMAL CURRENT FORMULA AND ONSAGER RELATION 

Transport through a QD coupled to two reservoirs in the presence of external voltages and temperature gradient 
between two reservoirs can be described by the Anderson single impurity model: 

^=Y1 ^nkAka^nka + Erf X! 4aCda + Uud^Tldl + ^ (^r/^fe^C^^ + H.c), (1) 

where e^feo- represents the conduction electron energy of the lead r], and is the discrete energy level in the QD. 
'^■qka (cr/fecr) are the creation (annihilation) operators for electrons in the lead "q {— L, R), while c|^^ {('da) for electrons 
in the QD. When a total external voltage V is applied between the two leads, their chemical potential diffi'erence 
is — A^R = eV. The two leads are assumed to be in local equilibrium with respective temperature and their 
distribution functions are given by f-qiui) = [1 + cxp (o) — p^ri) / k^T.q\^^ {rj =L or R). This assumption is practically 
correct because the two reservoirs respond to an external perturbation much faster than the center region, i.e., 
the QD. The other parameters U and V^, stand for the Coulomb interaction, and the coupling between the QD 
and the reservoir rj, respectively. For the transport problem concerned in this paper, electrons in the QD are in 
a nonequilibrium state, to be determined by their coupling to two leads and by the applied voltage. In order to 
describe the nonequilibrium state of electrons, we define the retarded (advanced) and lesser (greater) GFs for the QD 
as follows: G:,i^\t,t') ^ ±^0{±tTmc,^{t),cl{f')}), Gt{t,t') ^ ^{cl{t')c,^{t)) and G>,(t,t') ^ -*(c,,(t)4(t')>- 
The particle current and energy fiux Jsrj fiowing from the lead rj to the QD can be evaluated, respectively, from 
the rate of change of the electron number operator Nj^{t) = J2k a '^lika(^)%ka^^) ^^^^ °^ change of energy 

operator Hj^{t) = 'Ylik^a ^rika^ljkcr ityr,kai^) of tlie lead 77:^^ 

k.(T k,(T 

^^"^'^ = = Aa^'^'^''^''^'^^^*^''''^) = 'i{12'vk.KikAt)CdAt) - Kr4ws,.(o]), (3) 

k,<7 k,<7 

which involve the time-diagonal parts of the correlation functions: G'^cr r/ka^^'^') = ^i^rikcr(^')'^dcri^)) ^^'-^'-^rikcrdcr^^'^') = 
^{^d<T^^')'^-nka^^)) ■ According to Ref. 24, the thermal flux Jq^ flowing from the lead 77 to the QD is determined as 

Jqij = Jeti — IJ'TjJtj- (4) 

With the help of the Langreth analytic continuation rules,^^ we obtain the following expressions: 
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(5) 
(6) 



in terms of the QD's GFs in the Fourier space. Here r^(a)) = 27r J2k a \^v\^S{^~^vkcr) denotes the strength of coupHng 
between the QD levc;! and the lead r]. 

In the presence of both a strong Coulomb interaction in QD and tunnelings between QD and leads, it is difficult 
to evaluate the retarded GF G^^ of QD accurately. Several approximation schemes have been proposed in literature 
to derive G^^, such as EOM with the decouphng approximation, SOPT, etc, where a retarded (advanced) self-energy 

for the interacting QD is written as S''^"-' = Sq''"'' + S^^"-*, with Sg^"'' = r^(w)/2 being the nonintcracting part 

coming from the tunneling of electrons from the impurity state to outside leads, and S^'"-* being the interacting part 

derived within these approximation approaches. The retarded (advanced) GF G^J'^^ thus has the form 
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Unfortunately, one can not straightforwardly get the lesser GF for the strongly correlated systems under out-of- 

equilibrium circumstance, as docs for the retarded GF. Several years ago, Ng proposed a simple scheme to obtain the 
lesser GF G^^ from the retarded and advanced terms in order to study ac Kondo resonances in nonlinear transport 
basing on EOM approach. He assumed that S< (uj) = AT^^ (lu) and (w) = AEq (w) where A is an unknown function, 
while Sq = i'l2ri^ni^)fni^) ^^'^ ~ ^'^'l2rj^vi^)\-^ ^ fvi^)] ^^'^ nonintcracting lesser and greater self-energies. 
These lesser and greater self-energies are requested to satisfy the Keldysh relation S< — S> = S'' — S", leading to 
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or in explicit form, 



(8) 



(9) 



The lesser GF G^^ = E<|G^g.p is thus obtained. This is the central result of Ng's scheme. It has three advantages, 
initially addressed by Ng. that (i) it is exact in the equilibrium limit fi^ = M_r, (ii) it is exact in the nonintcracting (U = 
0) limit under general nonequilibrium situations, and (iii) the continuity equation Jhit) = —Jnit) is automatically 
satisfied in the steady state limit. ^ With the help of this ansatz, as long as one obtains the retarded GF properly 
describing the strongly correlated system from certain approximative method, the lesser GF can be derived, and thus 
the transport problem can be investigated. By means of Eq. (9), we easily obtain 



Jr, = -\j d^n^) [fni^) - ^(^)] ImGS<,(a;), 
Jqt, = ~\j dujT{w){uj - 11^) [fr,{uj) - /^(w)] ImGS„(t^), 



(10) 
(11) 



where r(a;) = T l[ij)T ii(ijo) / \r l(u)) + Tii{ijj)] and fj ^ rj. Note that here we arrive at exactly the same current formula 
(10) as that in Ref. 25 without introducing the assumption of a proportional coupling Tl{lo) = \Tfj{uj) (A is a 
constant). Since the pioneer work of Ng,^ several attempts have been made to generalize this ansatz to study Kondo- 
type transport in more complicated devices containing interacting QD, such as normal metal-QD-superconductor^* 
and superconductor-QD-superconductor^^. Nevertheless, the validity of this ansatz is worth being further examined. 
As mentioned in the introduction, verification of the Onsagcr relation is a natural choice for this purpose. 

The Onsager relation is concerned with the linear response of the particle current [Eq. (10)] and the heat flux 
Jq^ [Eq. (11)] driven by small bias voltages — Hfj = Sfj, and small temperature gradients — T^ = ST: 
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with T{u) = T{u) ImGJJ^(a;)5^=o,5T=o and 



'11 ~ .... .... 



dLoTico) ( ) , C,2 = — I dLoTico) [ ) , (14) 

T 
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/i 7 ^ ' \ djj. ) h J ' ' \ dT 

£21 = / du;T(a;)(a; - m.) (^^)^ > = ^ / ^'^^('^)('^ " /^.) (^^) • (15) 

We can easily see that Ng's ansatz Eq. (9) preserves Onsager relation £12 = £21 automatically. Furthermore, the 
result is independent of the approximation adopted in deriving the retarded GF. 

III. THERMOELECTRIC EFFECTS IN THE PRESENCE OF KONDO CORRELATION 

As shown in Eq. (12), both the bias voltage and the temperature gradient between two reservoirs can give rise to 

2 

particle current. The current induced purely by a small bias voltage reflects the electric conductance G = —^Cn, 
while the thcrmopowor 5* measures the voltage difiference needed to eliminate the current due to the temperature 
gradient between the leads, given in linear regime by 

In this situation, we can simplify the thermal flux Eq. (13) as Jg^ = —k6T, in which 
is the thermal conductance. 

Comparing the explicit expressions of £11 and £12 [Eq. (14)], we can roughly address that in low temperatures, the 
electric conductance G is determined by the transmission probability, or the density-of-state (DOS) of QD, at the 
Fermi energy of the leads, while the linear thermopower S depends sensitively on the energy dependence of DOS, 
implying that it contains information different from the electric conductance. So far, the thermopower of QD in the 
kondo regime is still much less studied in literature. In this section, we attempt to numerically calculate the linear 
Kondo-type thermoelectric effects in QD. The remaining task is to choose a suitable approximative scheme, which 
can provide; the retarded self-energy [or DOS p{i^')] capable of properly describing the Kondo physics in a wide range 
of parameters of interacting QD. In the present paper, a modified SOFT developed in Ref. 8 is adopted, in which 

with a = n(l — n)/no(l — no) and 6 = (1 — 2n)/no(l — no)U. n is the occupation number of the QD level which 
should be determined self-consistently. E^^^ and no are, respectively, the second order self-energy in U and a fictitious 
particle number, both of which are obtained from the bare GF Gq = — ed — Un — Eq). 

In actual calculation, an assumption that the tunneling strength is independent of incident energy is taken in the 
wide-band limit and a symmetric system Vl{w) = Ti{{u)) = F is focused. In the following we take the coupling 
strength F as the energy unit and the Fermi level of the lead to be zero. 

In Figs. 1(a)- (d) we plot the equihbrium DOS p{oj) for the QD with U = 7 and several different energy levels 
Cd = —2, —3.5, —5 and 2 at different temperatures. We can observe that the Kondo resonance peak in DOS is clearly 
resolved for the Kondo systems ed = —2, —3.5 and —5 at low temperature T = 0.01, but nearly dispeares at high 
temperature T — I. Of course there is no Kondo peak in DOS for the non-Kondo system ed — 2. In Fig. 2, the 
calculated electric conductance G (a), thermopower S (b) and thermal conductance n (c) are displayed as functions 
of the gate voltage, i.e., the energy level in QD. As expected, the conductance demonstrates a single peak structure 
and nearly unitary limit near symmetric point Cd = —U/2 at very low temperature, splitting peaks with increasing 
temperature and a minimum at = —U/2 for very high temperature. It manifests the equilibrium DOS of QD at 
the Fermi energy of the lead, p(0), at low temperature, but not its concrete shape of the Kondo peak: symmetric or 
nonsymmetric (right or left of center) around the Fermi energy of the lead. In order to detect DOS in the whole range 
of energy, one has made an attempt to measure the differential conductance of QD. However, the finite bias voltages 
between two leads result in a large change of DOS, giving rise to a splitting of the Kondo peak. As mentioned above, 
an alternative way of addressing this problem is to explore the thermopower, because it is relevant with the product 
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of the incident electron energy lo and the DOS of QD p{lo). For example, the perfectly symmetric shape of DOS 
aroimd a; = at the symmetric case td = —3.5, as shown in Fig. 1(b), results in exactly zero of the thermopower; 
while the positive (negative) thermopower is attributed to slightly slanting of DOS towards the left (right) for = —2 
(—5) [Figs. 1(a) and (c)], implying an oscillatory behavior around zero which can be controlled by the gate voltage. 
We also observe from Fig. 2(b) that the magnitude of oscillating in thermopower is largely decreased with lowering 
temperature due to the Kondo-suppressed deviation of DOS from the symmetric shape. Fig. 2(c) reveals that the 
thermal conductance k has similar behavior as the electric conductance G at low temperature, but quite different at 
high temperature T = 1. 

Fig. 3 displays the temperature dependence of these three quantities. We observe from Fig. 3(b) that the ther- 
mopower shows a logarithmic riseup with increasing temperature, leading to a broad maximum, and subsequently 
decreases in the high temperature regime. But temperature increasing can not cause a change of the sign of the 
thermopower, except for very high temperature where the Kondo effect is quenched. Fig. 3(c) shows that the thermal 
conductance has similar temperature dependence for these different systems. As a result, the thermal conductance is 
not a suitable tool to explore the Kondo effect. The interesting physical quantity is the ratio between the thermal 
and electric conductance. The classical theory yields that thermal and electric transport in bulk metals satisfy the 
Wiedemann- Franz law k/TG = 7r^/3e^. In mesoscopic system transport occurs through a small confined region, and 
consequently does not satisfy the Wiedemann-Franz law in general. We depict this ratio in Fig. 4. Surprising recovery 
of the Wiedemann-Franz law is observed at very low temperature, where transport through QD is dominated by Kondo 
correlation, implying that the Landau Fermi liquid state is rebuilt in this situation. But substantial deviation from 
the classical value emerges with increasing temperature due to disappearance of the Kondo effect. This temperature 
behavior provides a good trademark for exploring the onset of Kondo correlation. 



IV. CONCLUSION 



We have studied the particle current and thermal flux through interacting QD on the basis of the nonequilibrium 
Green's function approach and Ng's ansatz. The advantage of the ansatz is its capability of evaluating the lesser 
(greater) GF from the retarded and advanced GFs derived in certain rational approximation scheme. The validity 
about the ansatz and the resulting linear transport coefficients have been examined in term of Onsager relation. We 
have also emphasized that the same electric current formula as in Ref. 25 can be obtained without the assumption of 
proportional coupling Tl{ui) = \Tji{ui). 

In the wide band limit, a modified SOFT has been employed to calculate the retarded GF for interacting QD and 
with this retarded GF the thermoelectric effect in the Kondo regime has been investigated. We have found that 
the thermopower exhibits a oscillatory behavior around zero due to nonsymmetric shape of the Kondo peak in the 
DOS, giving rise to a change of sign of the thermopower which is controllable by tuning the gate voltage. These 
results demonstrate that measuring thermopower can provide useful information of the DOS in QD. Furthermore, 
our calculation reveals that the magnitude of oscillation is largely suppressed by the Kondo effect. Finally, we have 
explored the temperature characteristic of the thermoelectric effect and predicted that at low temperature regime 
thermal transport satisfies the classical Wiedemann-Franz law, which can be taken as a trace of the Kondo correlation. 
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FIGURE CAPTIONS 

Fig. 1 The DOS of the interacting QD U = 7 calculated from the modified SOPT with several energy levels ea = —2 
(a), —3.5 (b), —5 (c), 2 (d) and different temperatures T = 0.01, 0.5, and 1. In this and following figures, F is 

chosen as energy unit. 

Fig. 2 (a) The electric conductance G, (b) the thermopower S, and (c) the thermal conductance k, as functions of the 
energy level for several different temperatures T = 0.01, 0.5, and 1. The on-site Coulomb interaction in the QD 



is ?7 = 7. 



Fig. 3 (a) The electric conductance, (b) the thermopower, and (c) the thermal conductance vs temperature for the QD 
with U = 7 and = —1, 0, and —5. 

Fig. 4 The ratio factor between the electric conductance and thermal conductance vs temperature for the QD with the 
same parameters as in Fig. 3. 
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